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Abstract 

kjr)! We develop a method to compute the sphaleron rate in the electroweak broken phase non- 

perturbatively. The rate is somewhat slower than a perturbative estimate. In SU(2)xU(l) 
Higgs theory at the physical value of 0^, and assuming that the latent heat of the phase 
. transition reheats the universe to the equilibrium temperature, baryon number erasure after 
^ ! the phase transition is prevented only when (X/g 2 ) < 0.036. 
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Cl,; 1 Introduction 

D \ In the past ten years the effort to understand electroweak baryogenesis has driven the study 
of the physics of the hot electroweak plasma at and around the electroweak phase transition. 
The relevant physics has been put on a much sounder foundation, and in particular we now 
^ ■ possess reliable nonperturbative studies of the strength of the phase transition and of its 
disappearance at sufficiently high scalar (Higgs) self-coupling 0, 

Electroweak baryogenesis in the minimal standard model is now ruled out on two inde- 
pendent grounds. First, it cannot produce enough baryons because there is not enough CP 
violation in the minimal standard model §]. Second, even if it could produce the baryons, 
they would be wiped out in the aftermath of the phase transition, because the rate of baryon 
number violating processes does not shut off fast enough after the phase transition. This is 
because it is excluded that the standard model Higgs boson can have a weak enough self- 
coupling to supply the required strength to the phase transition Q. In fact, experimental 
bounds on the Higgs mass are now high enough that if the minimal standard model is the 
right physical theory, we can conclude that there was no cosmological electroweak phase 
transition. 

Extensions of the standard model, such as the minimal supersymmetric standard model 
with a light stop, appear to be alive and kicking for at least part of their parameter space. 

1 e-mail: guymoore@physics.mcgill.ca 
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Figure 1: "Cartoon" of the free energy dependence on Nqs- 



To determine what parts of their parameter space are viable, we demand that they pass 
two tests. There must be enough CP violation to produce sufficient baryon number to ex- 
plain astrophysical observations; and baryon number must not be subsequently erased. The 
production depends on where CP violation appears in the theory and on the mechanism by 
which it causes a baryon number surplus. We will not address this question here. The erasure 
depends on the efficiency of baryon number violating processes after the phase transition. 
This paper will study this erasure nonperturbatively. 

Arnold and McLerran established the basic picture for baryon number violation in the 
electroweak broken phase 0. Baryon number change in the electroweak theory is propor- 
tional to change in the Chern-Simons number (Nqs) of the SU(2) fields, because of the 
anomaly ||. But in the broken phase, there is a barrier to changing A*cs- To change from 
one integer value to another, one must pass through a half integer state, and while there are 
configurations of integer iVcs with arbitrarily small energy, the energy of a configuration of 
half integer Nqs is bounded by the energy of Klinkhamer and Manton's "sphaleron" con- 
figuration 0. At finite temperature the same statements are true for free energy, and one 
can use Langer's formalism || || to find the rate of baryon number violation. The classic 
cartoon of free energy as a function of Nqs is shown in Fig. |TJ. The hard part of getting 
from, say, Nqs = to Nqs = 1 is getting over the peak (the sphaleron barrier) in between. 
Let us assume that, whenever the system passes over the peak, it subsequently gets stuck 
in the new minimum. Then the diffusion constant for Nqs is set by the flux of states in the 
thermal ensemble going over the barrier. The diffusion constant for Nqs, 



is the probability per unit time of crossing the peak, which is the flux of states in the thermal 
ensemble over the peak. The flux of states over the peak is the probability to have iVcs within 
e/2 of the peak, times the mean value of \dNcs/dt\ measured at the peak, divided by e. Since 
the diffusion constant should be extensive, we generally refer to the diffusion constant per 
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unit volume, 1^ = 7d/V, often called the sphaleron rate. 

The height of the sphaleron barrier in the immediate aftermath of the electroweak phase 
transition depends on the broken phase expectation value of the Higgs condensate, <fi = 
\J (0 2 ), at the ambient temperature. In the minimal standard model the latent heat of the 
transition is not quite enough to heat the plasma back to the equilibrium temperature, but 
in super symmetric extensions, it may be; so we need the jump in ((f) 2 ) at T c . We now have 
reliable nonperturbative studies of this quantity JT|, [|. These studies were motivated largely 
so the determined value of (0 2 ) could be plugged into the calculation of the sphaleron rate. 
However, the present state of the art in the sphaleron rate is a one loop calculation. At 
the equilibrium temperature, we know that perturbation theory for infrared quantities such 
as is at best an expansion in X/g 2 , and we know of no reason why the expansion for 
should be better behaved than the expansion for (4> 2 ). At tree level, the phase transition is 
second order, and T d is large just below T c 0. At one loop, the phase transition becomes 
first order and is exponentially suppressed up to T c . To see how well perturbation theory 
converges, we would like a two loop calculation. Our experience from the calculation of (4> 2 ) 
is that two loop corrections may not be small. 

Unfortunately, extending the perturbative calculation of Td to two loops is immensely 
harder than it is for (<p 2 ). Both are background field calculations, but the background field 
in the effective potential calculation is spatially uniform, and the diagrams can be computed 
using Fourier techniques. The sphaleron background, in contrast, lacks translation symmetry, 
and is only known numerically. The one loop calculation demands summing the eigenvalues 
of numerically determined fluctuation eigenmodes in this background. At two loops one 
needs to compute mutual interactions between the fluctuations, which will involve double 
and triple sums over the eigenmodes of overlap integrals; and the interaction Hamiltonian 
depends on the numerically determined sphaleron background. In our opinion it is actually 
easier to make a direct nonperturbative calculation of than to perform this two loop 
calculation. 

In the symmetric phase we can now determine with fair reliability, using real time 
techniques [11], [L2JQ, but this does not work in the broken phase, because the rate is 
too small for a transition to occur in a reasonable amount of Hamiltonian time. But if we 
knew what we meant by the horizontal axis in Figure |l|, then we could use nonperturbative 
lattice tools to find the diffusion constant. Both the probability to be within e/2 of the 
peak, and the instantaneous mean value of \dNcs/dt\, are thermodynamic properties, and 
while the probability to be near the peak is exponentially small, we could measure it using 
multicanonical Monte-Carlo techniques. That would give up to an order unity dynamical 
prefactor, reflecting the fact that when the system crosses the peak, it does not necessarily get 
stuck in the new phase, but may turn around and re-cross immediately. In the broken phase 
we expect this dynamical prefactor to be on order 1; we will estimate it in the conclusion. 

Our approach in this paper is precisely to find a way to make sense of Fig. [I] in a nonper- 
turbative lattice study of the dimensionally reduced electroweak theory, and to construct the 

2 Recently it has been pointed out jlj], correctly, that in the parametric limit of small a w there are loga- 
rithmic corrections to the functional form assumed for fitting in [fl2| . The numerical value of the coefficient 
of the logcan be computed in a simple effective theory fll4[] , and the result is that the coefficient of the log 
is small |lq| ; so while there is an uncontrolled systematic in the most recent results for the symmetric phase 
rate JT^] , the error caused is relatively small. 
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figure and find the probability to be at the peak by multicanonical Monte-Carlo techniques. 
We also need to measure (\dN cs /dt\) at the peak, meaning the value averaged over the en- 
semble restricted to configurations with Ncs = 1/2. This paper will present the crux of the 
technique, with the essential ideas; but we will postpone some technical details to a longer 
sequel . The outline of the paper is as follows. First we construct a suitable definition 
of Ncs- Then we use it to perform a multicanonical Monte-Carlo determination of IV We 
compare our results at three values of x = X/ g 2 to results of a perturbative calculation of 
Td made using the two loop effective potential. We determine that, in the minimal standard 
model, or any theory where heavy degrees of freedom can be integrated out until it looks like 
the minimal standard model |TJ|, the sphaleron bound at the physical value of the Weinberg 
angle is x = \/g 2 ~ 0.036. 



2 Definition of Nqs 

We will treat the thermodynamics of the hot standard model in the dimensional reduction 
approximation, so the configuration space and probability distribution are the same as ther- 
mal, classical Yang-Mills theory with certain Higgs mass counterterms. We use notation 
and field normalizations appropriate for the 3+1 dimensional classical field theory and will 
determine the instantaneous value of \dNcs/dt\ in the classical theory context. But we need 
an appropriate definition of Ncs- 

Two features are essential to the definition of Ncs'- 

1. iVcs is the integral of the total derivative (g 2 /%tt 2 )E^B^; and 

2. iVcs is an integer for a vacuum configuration, Bf = 0. 

To determine the rate T d by constructing Fig. [I] we only need to define N cs modulo 1; so 
we can set Ncs = for all vacua, and we may hope for a gauge invariant definition. 

To measure Ncs of some three dimensional gauge field configuration we then find a path 
from that configuration to a vacuum and integrate EfBf along that path. The most natural 
path to choose is the cooling path; following |HJ we define a cooling or dissipative time r, 



under which gauge fields evolve by dissipative dynamics 
dA1{x,r) dH YM (A(r)) 



dr dA?(x,r 



A?(x,0)=A«(x). (2) 



We write in continuum notation for clarity; the covariant lattice implementation is in [fLl 
and Hetrick and de Forcrand have used a roughly equivalent procedure to solve the Gribov 
gauge fixing problem ||17|| . We cool only the gauge fields, under the Yang-Mills Hamiltonian, 
since Ncs should be a function of gauge fields alone; not cooling the Higgs fields also avoids 
some technical problems. 

In a finite volume, the gauge fields will cool towards a vacuum off a measure zero subspace 
of configurations which cool towards saddlepoints. Hence we can define Ncs as 
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N cs = £^l <h j<f\vE';(x.T)B'!{.v.T). (3) 
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where E£(t) is the field strength in the i,r direction. The standard lattice implementation 
of EfBf |T3 is not a total derivative ||19||, but it is to good approximation for smooth fields, 



by which we mean fields with almost no UV excitation. Cooling removes UV excitations 
much more efficiently than IR ones-in a linear theory the amplitude of a mode decays as 
exp(— rk 2 )-so the shortcomings of the definition of EfBf, and regulation sensitivity more 
generally, go away at a very small value of r, and are not a problem. Also note that the 
cooling procedure and the definition of E^B? are gauge invariant, and so is our definition of 
N C s- 

We cut off the small r part of the integration to remove UV contributions to Nqs which 
are unrelated to winding number change. In the continuum abelian theory, where 

2 

iVcs = Jr^ [ d z xe ljk F ij A k , (4) 
the mean square value of Nq§ we would get if we used To = is 
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(A^cs) = Y^f^I J ^xd^ye^e^F^A^F^A^y)) . (5) 

Using Wick's theorem and the momentum representation of the propagator, in Feynman 
gauge, this becomes 

(Age) = T7^ f***V /^^ +gH ^^n { ^ T ]Y kn + 
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so Nqs will be Gaussian distributed with a linearly divergent variance. On the lattice, the 
UV divergence will be cut off by the lattice scale; the coefficient was found by Amjorn and 
Krasnitz p| and is (A^ s ) = (1.44 x 10" 5 )^ 4 UT 2 /a. 

SU(2) theory will have similar UV contributions, since in the UV it looks like three copies 
of the abelian theory. This UV contribution has nothing do with topology, since it appears 
already in the abelian theory. If we used r = then the measured value of iVcs would be 
an "interesting" IR piece plus this nontopological UV piece. The probability distribution of 
A^cs would be the convolution of the probability distributions for the two, and if the noise 
distribution were broad, this would spread Nq$ ~ configurations out to dominate the 
sample at iVcs = 1/2, destroying the peak in Fig. |l| Only if the noise distribution is narrow 
will the probability distribution be undistorted. 

At finite tq, and including the mass for the gauge fields which appears in the broken 
phase, the UV contribution to iVcs, Eq. (^), becomes approximately 

c/ 4 T 2 U r d 3 p p 2 e - 4 P 2T0 



^CS) " J (2^3 (^2 + m 2 w) 2 ■ ( 7 ) 

(Note that r has units of length squared, not length.) Choosing tq > l/m 2 ^ controls the noise 
sufficiently, while leaving Nq$ with the feature which is essential to our endeavor, namely 
that a configuration which is very nearly a sphaleron (meaning that it cools to very close to 
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the "sphaleron" saddlepoint, but then slips off and finds a vacuum) will have iVcs very close 
to ±1/2. We have checked that the part of the cooling path from (almost) the sphaleron 
saddlepoint to the vacuum gives (g 2 /8% 2 ) J E^Bfd 3 xdT of almost exactly ±1/2. 

This definition of Ncs has a somewhat arbitrary parameter To, so we cannot give a 
direct physical interpretation to the probability distribution of Ncs- Also, the probability 
distribution of Ncs will depend on r . Consider a configuration which is "fairly close to" 
a sphaleron but which is already rolling away from the saddlepoint at cooling time r . If 
we increase r , it will roll further out of the saddlepoint before we start to integrate E ■ B, 
and if we decrease To then it will be closer to the saddlepoint. Hence, choosing a larger 
To will spread out the determined values of Ncs of similar, nearly sphaleron configurations 
more widely, diluting the probability distribution near Ncs = 1/2 (and concentrating the 
probability distribution near Ncs = 0). This reduces the population of configurations with 
iVcs within e of 1/2. However, spreading out the configurations increases dN cs /dt, evaluated 
on an ensemble of configurations with N = 1/2, by "almost" the same factor, and so the rate 

"almost" does not depend on tq. By "almost" we mean that this statement is strictly true 
when we can neglect the contribution to \dN/dt\ arising from residual UV fluctuations. The 
more we cool, the larger the contribution to dN/dt from motion in the unstable direction 
and the smaller the contributions from the stable directions, so the truer this becomes. In 
a complete calculation including the "dynamical prefactor" this remaining tq dependence is 



also accounted for For this work we will use a value of To where we can directly check 
that fluctuations about the sphaleron are largely absent for configurations with Ncs = 1/2. 

We make two further approximations to make measuring Ncs numerically practical. The 
cooling is very numerically expensive, but after cooling by r > a 2 all UV information on 
the lattice has been destroyed, and cooling on such a fine lattice is redundant. We set up a 
coarsened lattice of all even sites and continue the cooling on it. We have compared the value 
of Ncs obtained in this way to the value without coarsening, and aside from a slight and 
measurable renormalization, the difference is a tiny amount of noise. We also extrapolate 
the large r exponential tail of / E ■ Bdr to avoid having to run r to arbitrarily large times. 



3 Multicanonical measurement 

Now that we have a definition of Ncs (modulo 1), we want to measure the probability to 
have a particular value, that is 

P(x) = -J- / VVDAeT f)H (*AQ{Nc& - x)Q(x + 5x - N cs ) , (8) 
Zox J 

where Z is the same integral without the O functions, which select out the range of iVcs 
between x and x ± Sx. We are making the dimensional reduction approximation including 
the integration over the A Q field, though we do include the U(l) hypercharge field at the 
physical value of Qw the Weinberg angle, using the noncompact lattice implementation. 

To determine P(x) accurately even where it is small, we perform a multicanonical Monte- 
Carlo calculation. That is, we rewrite Eq. (j8|) as 

P(x) = Ybi\ (v<S>VAe-P H ^^e fiNcs) ) e~^ Ncs) Q(N cs - x)Q(x + 8x - N cs ) , (9) 
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Figure 2: Free energy (left) and (^broken — ^sjmm) / (9 2 ^ 2 ) as functions of N cs at x = X/g 2 = 
0.039, in a (lQ/g 2 T) 3 volume. The plot of (ft 2 is a check that the volume used was large 
enough to prevent the sphaleron from causing a transition to the symmetric phase. 



where f(x) is a function chosen by bootstrapping to be approximately equal to — lnP(x). 
We do the integral by finding a Markovian process which samples configuration space with 
weight e~ )3H ^' A ^e^ Ncs ^'D^T)A, and replacing the integral with a sum over a sample generated 
by that Markovian process. Since Nqs is a complicated function of the configuration, this is 
easiest done by finding a Markovian process which samples with canonical weight and adding 
a Metropolis accept reject step to account for eA 

We generate the canonical Markovian process with the algorithm from [|l[], extending it 
to include the noncompact U(l) subgroup, and reducing the efficiency of the gauge field 
updates, to bring the Metropolis accept rate up to about 50%. To measure \dNcs/dt\, we 
augment the system with Gaussian "momentum" degrees of freedom, pick momenta from the 
thermal distribution, and evolve all fields forward under Hamilton's equations (implemented 
by leapfrog) for a (very short) time. We determine dNcs/dt from the difference between 
iVcs of the initial and final configurations. This differs from a Hamiltonian evolution of the 
classical system because we do not enforce Gauss' law; but since B is transverse and Gauss' 
law only affects the longitudinal E field, the instantaneous value of E ■ B should not be 
affected. 

We have measured T d for three values of scalar self-coupling, x = (X/g 2 ) = 0.047, 0.039, 
and 0.033. We used a 40 3 box with periodic boundary conditions and a lattice spacing of 
a = 2/(hg 2 T) (f3 L = 10) for the larger values of x and a = l/(3g 2 T) (/3 L = 12) for x = 0.033. 
We determine the equilibrium temperature as in f20fl . The lattice action is 0(a) improved 
using the relations in [21, |22j , so finite lattice spacing systematic errors will be smaller than 



statistical errors. The very large volume is necessary to prevent the tails of the sphaleron 
from meeting because of the periodic boundary conditions. It is also necessary because the 
Higgs field condensate has a zero at the core of a sphaleron, so the sphaleron looks somewhat 
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860 


920 


1000 


nonperturbative 


0(T c )/^T c 


1.38 ±0.02 


1.60 ±0.01 


1.82 ±0.03 




- ln(r d T" 4 ) 


24.7 ±0.4 


28.3 ±0.4 


31.2 ±0.6 


mixed 


- ln(r d T" 4 ) 


22.6 ±0.3 


25.9 ±0.2 


29.3 ±0.5 



Table 1: Perturbation theory versus nonperturbative IV Appearances of T c ~ 4 are really 
(2.5g 2 T c ) -4 , and d/dy means g^d/ 'd{m\/T 2 ). The error bars for the nonperturbative 4>(T C ) 
are dominated by statistical errors in the determined value of T c ; errors in the nonperturba- 
tive value of T are statistical errors from the Monte-Carlo. The "mixed" results use the two 
loop value of B but the nonperturbative value of <f)(T c ). 



like a symmetric phase bubble; the large volume prevents the sphaleron from stimulating a 
transition to the symmetric phase. 

We plot the numerically determined free energy and (0 2 ) versus Nqs for the x = 0.039 
data in Fig. g. The value of (\dN CS /dt\) at N cs ~ 0.5 was (0.22±0.05)a- 1 = (0.55±0.13)# 2 T. 
Recall that both this value, and free energy distribution, depend somewhat on our choice of 
Tq, which was r = 22.5a 2 = 3.6/g A T 2 . With less cooling, the free energy would rise more 
evenly and would be less flat on top; so the shape of the free energy curve should not be 
overinterpreted. The determined rate T d should be independent of r , though. 

We present the results in Table |l|. They should be compared to the symmetric phase 
rate, which, using g 2 = 0.40, is T d = (29±6)a^T 4 = exp(-13.9±0.2)T 4 As expected, 

the rate is orders of magnitude smaller in the broken phase than in the symmetric phase. 



4 Comparison to perturbation theory, erasure bound 

Next we compare the rate to a perturbative estimate. One loop perturbation theory gives 

r, = (io) 

Here <fi is the broken phase Higgs condensate expectation value, uj- is the unstable frequency 
of the sphaleron, Mt r MV ro t are zero mode factors, k is the one loop fluctuation determinant, 
and -E sp h is the energy of the Klinkhamer Manton sphaleron, using the tree level Hamiltonian. 
For small scalar self-coupling, — T In k equals the energy due to the one loop effective potential 
term, plus a modest correction [f25|. We can guess that the dominant two loop corrections to 
Eq. (|10D are absorbed by including the two loop effective potential terms in the Hamiltonian. 

3 Again, this result misses a logarithmic correction which is computable but not computed 
4 The definition of T used in || ^3) is the response rate to a chemical potential, which is half the diffusion 
rate ||; so Eq. © differs by a factor of 2 from the expressions in those references. 
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So it seems reasonable to estimate the sphaleron rate by Eq. (0), but setting k = 1 and 
solving for the sphaleron energy using the two loop effective potential at the equilibrium 
temperature. One should also solve for the zero modes and u;_/</> at this value, but they 
are very weak functions of the effective potential [p^| . We use the values from [^3| at 
x = {X/g 2 ) = 0.04 for these, but solve for the sphaleron energy, i? sp h = 4nB(j)/g, numerically, 
using the two loop effective potential at T c . We use the two loop potential presented in [p6| . 
without pieces from longitudinal gauge bosons (assumed integrated out). We also drop two 
loop terms proportional to Xg 2 or A 2 , because the perturbative determination of <fi is an 
expansion in X/g 2 , and such terms contribute at the same or higher order as unknown 3 loop 
terms. (Including those terms moves <fi closer to the nonperturbative value.) We compare the 
results to the numerically determined values in Table |I[ The "two loop" analytic sphaleron 
rate is about exp(2.5) times faster than the numerically determined nonperturbative rate. 
The difference is more than can be explained by the difference in 0. 

We should compare the sphaleron rate to the limit set by requiring that baryon number 
not be erased. The rate at which sphalerons degrade baryon number is [f| 

-^r,T-3, (n) 



iV R dt 



where Np = 3 is the number of families, and the numerical factor 13Np/4 would be smaller 
in theories, such as supersymmetry, in which additional degrees of freedom can store baryon 
number.^ Integrating from the end of the phase transition to the present day, 

1 3 Nt? r°° 

In(JV B /JV B (r c )) = -W T d (T(t))T-\t)dt, (12) 

4 Jto 

where we have shown the dependence of on T and of T on t. 



Now lnT^ is very sensitive to J (0 2 ), and hence to T; so we can approximate lnT^T) 
lnT rf (T c ) + (T — T c )(d In T d /dT) \ t=t c , and perform the integral: 



ln(-M^B/JV B (r )))=ln(39/4)+lnf^ > l -J dh > T *FW 



T 4 ) \ Tdt 



T=T C 



By the chain rule, 

d In Yd d In T d dy dliaT 

Tdt dy dT dt 

where y = m^(T) / (g 4 T 2 ) is the dimensionless thermal Higgs mass squared. 
Two of these are easy. We get dy/dT from the 1 loop correction to [BC 

dy _ 8A + 4# 2 + # 2 (3 + tan 2 Q w ) 
dT ~ ' 

and we get d\n.T/dt from the Friedmann equation in a radiation dominated universe, 



(13) 



(14) 



(15) 



J_ = H 2 = 87rG7l2g * T i =► = - /^ull (iq) 

At 2 3 30 dt V 45 m p i ' 1 ' 



5 Again, there is a factor of 2 difference from the reference because they write in terms of the response to 
a chemical potential, which is half the diffusion constant. 
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where g* is the number of radiative degees of freedom in the universe (g* = 106.75 in 
the minimal standard model) and m p \ ~ 1.22 x 10 19 GeV is the Planck mass. Finally, we 
determine dlnT^/dy perturbatively, by varying y slightly from the equilibrium value and 
recomputing the two loop sphaleron rate. The dependence is quite strong. We include it in 
the table. 

The most widely cited discussion of baryon number erasure after the phase transition 
makes the approximation that the baryon number violation rate after the phase transition 
is constant for about one Hubble time fl2"7fl . In fact, because T d depends very strongly on y, 



which in turn depends strongly on T, most baryon number erasure occurs in the first 10~ 3 
Hubble times after the phase transition. Hence the initial rate of baryon number violation, 
Td{T c ), which prevents washout is 10 3 times larger than assumed in f2~7|| , leading to a weaker 
bound on Td(T c ), roughly 

- ln(r d (T c )T- 4 ) > 30.4 - ln(T c /100GeV) . (17) 

The values of g* and dy/dT are both larger in supersymmetric extensions of the standard 
model, by on order a factor of 2; so the bound, Eq. (|T7|), is weaker in those models by about 
1. Also note that, because Eq. (|T3|) is for the double log of Nb/Nb(T c ), failing to meet the 
bound by 1 means the baryon number is diminished by exp(exp(l)) ~ 15, and failure by 2 
reduces baryon number by exp(exp(2)) ~ 1600; so the bound is quite sharp. 

Interpolating between the values of x where we have measured, and including the estimate 
discussed in the conclusion for the dynamical prefactor, we get a bound of about x = 0.036 
in the standard model and x = 0.038 in the MSSM. 



5 Conclusion 

Our nonperturbative results show that the bound on x = X/g 2 is about 0.036 in the standard 
model. This result will be corrected only slightly in any theory where all non-standard 
model degrees of freedom are heavy enough to integrate out accurately. This should apply 
for instance to the minimal supersymmetric standard model unless the lightest stop is quite 



light | 16| . The bound weakens somewhat if the latent heat of the phase transition is not 
sufficient to reheat the universe to T c ; this should be checked in the MSSM. 

If additional degrees of freedom are light, then they can affect the thermodynamics of 
the sphaleron (aside from renormalizing couplings). The calculation should then be redone, 
including those degrees of freedom. However, we have shown that if a two loop perturbative 
calculation gets the value of <p about right, then its value for Yd is reasonably close. We do 
not expect light stop squarks to endanger this conclusion; while their contribution to the 
effective potential may be important, they only contribute to the scalar wave function at 
two loops, and they do not interact directly with the gauge fields at all. Hence, they should 
change the sphaleron bound mainly by changing the size of the Higgs condensate. We can 
paraphrase the results of this work as setting a bound on <fi of roughly <fi > 1.7gT, and this 
should be approximately the same in extensions with a light stop. 

Finally we should comment on the approximation that successive sphaleron crossings will 
be in uncorrelated directions. Arnold, Son, and Yaffe have recently criticized this assumption 
28fl , on the grounds that it ignores dynamical effects related to hard thermal loops. Their 
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picture has recently been borne out in the symmetric phase f!2l| . In the present context, we 
would expect the system to cross the sphaleron barrier multiple times before relaxing to a 
vacuum configuration if the plasma frequency u> p i is much greater than the unstable mode 
frequency uj-. It should be clear why; when uj v \ the electric field, and hence EfBf, 

will oscillate in sign, driving Nqs back and forth, on a time scale short compared to the time 
scale at which it "falls off" the sphaleron. On the other hand, if u p i < u;_, then we expect 
no such effect. A crude, parametrically correct estimate is that the rate we determined will 
be corrected by a factor of 



1 + 



\\g 2 T 2 " 



1 + 



i -i 



(0.47,? 



2 i2\-l 



Here we use the standard model value for uj 2 ^ 
0.03g 2 . Numerically, using g 2 



and the value of u 2 _ computed in |^3| for 
A ~ 0.03g 2 . Numerically, using g 2 ~ 0.4, this estimate reduces the rate by a factor of 2 for 
x = 0.033, changing the bound on x = X/g 2 by ~ 0.0015. The correction in the symmetric 
phase is larger, because the natural length scale of a winding number changing process is 
larger; if we estimate that scale as g 2 T/2, for instance, the correction would be ~ 7. It is 
possible to use the techniques of |12| to measure the prefactor numerically [fLlfl , but we will 
not do so here. 
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